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Abstract. Wc point out that the determination of the leptonic decay width of radially- 
excited quarkonia is strongly dependent on the position of the node typical of these 
>—} | excitations. We suggest that this feature could be related with the longstanding p~ it puzzle. 
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2: 1. INTRODUCTION 

i> : 

To study processes involving heavy quarkonia, such as decay and production 
mechanisms, in field theory, all the information needed can be parameterised by 
vertex functions, which describe the coupling of the bound state to its constituents 
and contain the information about the size of the bound state, the amplitude of 
probability for given quark configurations and the normalisation of the bound-state 
(U [ wave functions. 

^ ' In a parallel work 0, 0, HI], we have have considered production processes of 

J/ip, ip' and T at hadron colliders. We relied on a phenomenological approach with 
vertex functions, and after restoration of gauge invariance, we have reached an in- 
teresting agreement with production cross sections and polarisation measurements 
by CDF $ 0, 0| and PHENIX g . 

In the case of J/ip and T(1S), we have used the leptonic decay width to fix the 
normalisation of these vertex functions. This procedure has appeared as robust 
and reliable for IS states and is explained in the following. However in the case 
of radially excited states, an important feature has emerged from calculations: the 
decay width, and thus the normalisation, are strongly dependent on the position 
of the node appearing in the vertex function. This ambiguity in the determination 
of the decay width might also exist for hadronic decays, such as ip' pn, and is 
perhaps closely related to the p — ir puzzle. 
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2. OUR PHENOMENOLOGICAL APPROACH 



Our approach to build the quarkonium vertex functions does not rely on the Bethe- 
Salpeter equation [<| (BSE), usually used to constrain the properties of vertex 
functions; for example, to relate the mass of the bound state to the mass of its 
constituent quarks. Indeed, besides problems with gauge invariance, all predictions 
coming from BSE are realised within Euclidean space (see e.g. 

mM\m) and we 

have shown that the continuation to Minkowski space can be problematic 

We have therefore chosen to describe the quarkonia by a phenomenological 
vertex function and decided to remain in Minkowski space, at variance with other 
phenomenological models (see e.g. [3,I3|)- The price to pay for not using BSE is an 
additional uncertainty due to the functional dependence of the vertex. Fortunately, 
it can be cancelled if one fixes the normalisation from the study of the leptonic 
decay width, which is our main concern here. 



2.1. Choice for the vertex function 

It has been shown elsewhere (see e.g. [IT]) that other Dirac structures than 7 M 
for vector bound state are suppressed by one order of magnitude in the case of 
light mesons, and that this suppression is increased for the <fi meson. Therefore 
their effect is expected to be even more negligible in heavy quarkonia. 

As a consequence, the following phenomenological Ansatz for the vector meson 
vertex, inspired by spin-projection operators, is likely to be sufficient for our 
purposes. It reads: 

v,( P ,p) = r(p,p) 7/1J (i) 

with P the total momentum of the bound state, P = pi — p2, and p the relative 
one, p = (pi +p2)/2 as drawn on Fig. ^ This Ansatz amounts to multiplying the 
point vertex (corresponding to a structureless particle) by a function, T(p,P). 




FIGURE 1. Phenomenological vertex obtained by multiplying a point vertex, representing a 
structureless particle, by a vertex function (or form factors). 

The function T, which is called the 3-point vertex function, can be chosen in 
different ways. Two simple Ansatze are commonly used in the context of BSE. We 
consider both. They correspond to two extreme choices at large distances: a dipolar 
form which decreases gently with its argument, and a gaussian form: 

N p 2 
r (p,P) = ^ and T (p,P) = Ne^, ( 2 ) 

I 1 A 2 / 



both with a free size parameter A. As T(p,P) in principle depends on p and P, we 
may shift the variable and use p 2 — ^9 instead of p 2 which has the advantage of 
reducing to — \p\ 2 in the rest frame of the bound state. We refer to it as the vertex 
function with shifted argument. In the latter cases, we have -in that frame-, 

N -m 2 
r(p,P) = (i+|pT)2 and T(p,P) = Ne^. (3) 

2.2. Excited states 

As is well-known, the number of nodes in the wave function, in whatever space, 
increases with the principal quantum number n. This simple feature can be used 
to differentiate between 15 and 25 states. 

We thus have simply to determine the position of the node of the wave function 
in momentum space. To what concerns the vertex function, working in the meson 

rest frame, the node comes through a prefactor, 1 — , which multiplies the 

vertex function for the 15 state. Explicitly, T 2 s(p,P), for a node a node} reads 



V ^ode J (1 + |£) 2 V ^ode J 



e^~ 



(4) 



In order to determine the node position in momentum space, we can use two 
methods. The first is to fix a no d e from its known value in position space, e.g. from 
potential studies, and to Fourier-transform the vertex function. In the case of a 
gaussian form, this can be carried out analytically. 

The second method is to impose the following relation 1 between the 15 and 25 
vertex functions: 

\p\ 2 d\p\e=^ ( 1 - -$-) = 0. (5) 



^node 

The two methods give compatible results. 

3. NORMALISING: THE LEPTONIC DECAY WIDTH 

The width in terms of the decay amplitude Ai, is given by 



1 1 



, _ . 1 2 



2M (4tt 2 

where d2{PS) is the two-particle phase space [l5l |. 



M\ <k(PS), (6) 



inspired by the orthogonality between the 15* and 2S wave functions. 



The amplitude is obtained as usual through Feynman rules, for which we use our 
vertex function at the meson-quark-antiquark vertex. At leading order, the square 
of the amplitude is obtained from the cut-diagram drawn in Fig. El 
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FIGURE 2. Feynman diagram for 3 Si -> ££. 
In terms of the sub-amplitudes A^, B^ u and C pv defined in Fig- El we have 2 : 

J \M\ 2 d 2 (PS) = \^A^ B*£*jgL)C'*, (7) 

where the factor = (g^ — ) = J2i £ i,p £ i,p' resm ts from the sum over polar- 
isations of the meson and the factor ~ accounts for the averaging on these initial 
polarisations. 



3.1. Sub-amplitude calculation 

To what concerns the sub-amplitude B pp \ from the Feynman rules, and after 
integration on the two-particle phase space, we have 



BPP' =(ie) 2 



ttM 2 ( , pppp' 
nM 2 g p P '-8—(gPP +2 
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For A 11 ", from the Feynman rules and using the vertex functions discussed above, 
we have (see Fig. EJ) 



iA» u = -3e c 
= -3e c 



d 4 k 
d 4 k 



Tr (iT(k,P)Y 



i{^+\P + m) 
(k — y) 2 — m? + ie v ' ' (k + ^) 2 — m 2 + ie 



-(ie'y 



g^(M 2 + 4m 2 -4k 2 ) + 8k^k u -2P^P V 
( ' \{k-^) 2 -m 2 + ie)((k+^) 2 -m 2 + ie)' 



(9) 



cq is the heavy quark charge, —1 comes for the fermionic loop, 3 is the colour 
factor. 



In the Feynman gauge. The calculation can be shown to be gauge- invariant, though. 



k + ^P 




k-iP 



FIGURE 3. Feynman diagram for 3 Si — > 7*. 

We can easily motivate our choice of the argument of the vertex function, beside 
its simple relation with that of wave functions. Its main virtue is to regularise the 
integration in all spatial directions. In the meson rest frame, the k° integral can be 
done by standard residue techniques, and the remaining k integral is guaranteed to 
converge. If one uses k 2 instead of \k\ 2 as an argument, it is possible to show that 
along the light cone, one obtains logarithmic divergences, which would presumably 
need to be renormalised. Besides, the k° integral then also becomes dependent on 
the singularity structure of the vertex function. These two reasons make us prefer 
the shifted vertex in our calculation. 

A notable simplification can be obtained by guessing the tensorial form of A^ v . 
Indeed, current conservation (gauge invariance) for the photon can be expressed 
as: 



A^Py = 0. (10) 

It is equivalent to 

pp, pv 

iA^ = iF(g^ - -pp) = iFA^, (11) 

the coefficient F being -p since A^ = A^ u g v ^ = 3F. This assumption can be easily 
verified in the bound-state rest frame P — (M, 0,0,0) for which Eq. (|TU|) reduces 
to A^°M = A^° = 0. We shall check this at the end of this calculation. 

It is therefore sufficient to compute the following quantity, where we set \k\ 2 = K 2 
and define T(-K 2 ) = T(k,P), 



POO 

iA* = -3e Q / AnK 2 dKT(-K 2 ) x 
Jo 



dko 4(-2(fc 2 -^f)+4m 2 ) (12) 

^ (2nY((k-?r) 2 -m 2 + ie){{k+?r) 2 -m 2 + ie))' 



Let us first integrate on ko by residues. Defining E = v 1 K 2 + m 2 , we determine 
the position of the pole on k (still in the bound-state rest frame) as 

(k ± p) 2 - m 2 + is = (k ± - E 2 + is. (13) 



To calculate the integral on fc , we choose the contour as drawn in Fig. HJ Two 
poles —E — y and — E+ 4^ are located in the upper half-plane and the two others 
in the lower half-plane. We shall therefore have two residues to consider. The 
contribution of the contour Cr vanishes as R tends to oo. 
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FIGURE 4. Illustration of the contour chosen to integrate on ko- 

The procedure is correct only if the poles of the upper-plane are on the left side 
and those of the lower-plane are on the right 3 . 

nnintQ — T?. -1- , 

2 2 



This crossing would first occur between the points —E + 4f and E — M-. Solving 



for E, we have: 



M . / . / 2 , 1E v 

-E + — = E =4> m 2 = K 2 >0. (15) 

2 2 4 



If they cross, the integral acquires a discontinuity from the pinch and becomes complex. 



The crossing is therefore impossible when M < 2m. Therefore, to get a coherent 
description of all charmonia (resp. bottomonia) below the open charm (resp. 
beauty) threshold, we shall set the quark mass high enough to avoid crossing for 
all of these. This sets m c to 1.87 GeV {mo) and m;, to 5.28 GeV (me). Considering 
the variety of results obtained from potential models, this seems to be a sensible 
choice. 

Putting it all together, we get 

-3e Q f°° K 2 T(—K 2 ) (2K 2 + 3m 2 ^ 



71^ 



A* = — ^ / dK - v (16) 



VK^+m 2 (K 2 + m 2 -^) 



One is left with the integration on K for which we define the integral I depending 
on the vertex function whose normalisation is pulled off, 



dKK 2 T(-K 2 ) (2K 2 + 3m 2 
VK 2 + m 2 N (K 2 + m 2 - 



I^M,m)=l ,=== \, > ,l^Z, ' (17) 



J is a function of A through the vertex function T(—K 2 ) and is not in general 
computable analytically. In the following we shall leave it as is and express A^ u as: 

A"= ^-p-NI(A,M,m) =>A' W = ^NI(A,M,m)A^ . (18) 

TV 7T 

Finally, to what concerns the sub-amplitude C, we simply have 

<™V = (^V)t = ^V. (19) 

3.2. Results 

Now that all quantities in Eq. (J7J) are determined, we can combine them. Hence, 
we obtain 4 : 

J \M\ 2 d 2 {PS) = ^f=± (^NI(A,M,m)) 2 (w^lj A^A«V'A^' 

(20) 

The leptonic decay width eventually reads from Eq. ©: 

T "= N2 ^B^ M ^- (21) 



4 Recall that the projector A M „ satisfies A^A^ = 3 and A^A""' = A./ 



3.2.1. Numerical results 



The value of TV is in practice obtained by replacing T« by its measured value, eQ 
by y for c quark and ^ for b quark and, finally, by introducing the value obtained 
for / for the chosen value of A. We sketch in Fig. El and Fig. |U1 some plots of N for 
the J /if) and for the T(IS') to show its dependence on A and tuq and the value it 
actually takes. 
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FIGURE 5. 

(resp. left ). 



Normalisation for a dipolar (resp. gaussian) form for J /ip as a function of A: right 
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FIGURE 6. Normalisation for a dipolar (resp. gaussian) form for T(1S) as a function of A: 
right (resp. left ). 

TABLE 1. Set of values for m,, A and N obtained for 
the J/ip within BSE approach [Tq. 



m c (GeV) || Vertex functions 





Gaussian 


Dipole 


1.6 


V 


A (GeV) 


N 


A (GeV) 


1.59 


2.34 


1.33 


2.87 


1.7 


2.28 


2.08 


1.87 


2.54 


1.87 


3.31 


1.97 


2.76 


2.31 



We have chosen the range of size of the mesons, or A, following BSE studies |16j 
and other phenomenological models k| 3|> where the commonly accepted values 
for the J/ip are from around 1 to 2.4 GeV, and for the T(IS') from around 3.0 
to 6.5 GeV. For illustration, we put in Tab. Q several values of A obtained with 
BSE [HI, accompanied with the vertex- function normalisation. 



3.2.2. Result for i/j' and implication for the p-ir puzzle 

To what concerns the ip', we give the results for A=1.8 GeV, m c = 1.87 GeV 
and a gaussian vertex function (see Fig. EJ). The normalisation N' diverges for 
a no de — 1-35 GeV. This comes from the cancellation of the leptonic decay width for 
this value of a no d e - 

100 r M J 



c-10 



03 

CD 



0.1 1 

0.5 1 1.5 2 2.5 

a node ( GeV ) 

FIGURE 7. N' for the ip' as a function of a no d e - 

This cancellation can be traced back to the integrand of I (see Eq. (|17Jl ). which we 
note (see Fig. [HJ. The cancellation of the positive contribution by the negative 
one therefore occurs in this case for a no d e = 1-35 GeV. 

The constraint for the node position Eq. (jSJ) gives a no d e = 1-46 GeV for these 
values of A and m c . For this value of a no d e , the normalisation N' is 16.36. This was 
the value that we have retained for the numerical applications of 0. 

In the case of the hadronic decay ip' — ► pir, we expect it to proceed via an off- 
shell uj (see Fig. IHJ). The integral / involving the ip' vertex function will be slightly 
different and might be, for the same value of a no d e , drastically suppressed. 

We can therefore expect a severe suppression of the decay amplitude compared 
to the nodeless, i.e. J/ip, case. This suppression could in turn explains the p — n 
puzzle, namely that the measured ratio prfeiz^ is not of the order of 15 % 
(expected from the leptonic decays) but rather smaller than 1%. 
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FIGURE 8. Evolution of the integrand in / as a function of K and for three values of a no de 
in the tp' case. 




FIGURE 9. Possible decay of a tp' into a p - it pair via an off-shell u. 



4. CONCLUSION 

We have explained our approach to describe on a phenomenological basis the inter- 
nal dynamics of heavy quarkonia in the context of a Feynman-diagram calculation. 
We have also provided a way to represent the distinct features of 2S states, like 
the ip' . We have seen that a simple but robust constraint on the vertex functions 
for lS'-states could be achieved by imposing that their normalisation reproduces 
the leptonic decay width, through a simple leading order calculation. 

An interesting simplification can be also obtained by shifting the argument of the 
vertex function, namely the relative momentum of the quark inside the quarkonium, 
into a quantity that reduce to the tri-dimensional relative momentum in the meson 
rest frame. This enable us to work out analytically, for whatever vertex function, 
the integration on ko. 

Furthermore, we have shown that the leptonic decay width was very dependent 
on the node position, which, incidentally, is not a completely constrained param- 
eter. This induces the same effects on the normalisation (see Fig. EJ). This has 
consequences on production processes where our approach to describe internal dy- 
namics of heavy quarkonia can be applied (see l2,|3j]). It should be interesting 
to see whether this happens for other excited states. Finally, we suggest that this 
feature typical of radially excited states could be the awaited explanation for the 
longtstanding p — n puzzle. 
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